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Abstract 

We derive simple analytical formulae for the renormalization group running of 
neutrino masses, leptonic mixing angles and CP phases, which allow an easy un- 
derstanding of the running. Particularly for a small angle #13 the expressions 
become very compact, even when non- vanishing CP phases are present. Using 
these equations we investigate: (i) the influence of Dirac and Majorana phases on 
the evolution of all parameters, (ii) the implications of running neutrino parame- 
ters for leptogenesis, (iii) changes of the mass bounds from WMAP and neutrino- 
less double (3 decay experiments, relevant for high-energy mass models, (iv) the 
size of radiative corrections to #13 and #23 and implications for future precision 
measurements. 
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1 Introduction 



The Standard Model (SM) agrees very well with experiments and the only solid evidence 
for new physics consists in the observation of neutrino masses. Compared to quarks and 
charged leptons they are tiny, for which the see-saw mechanism provides an 

attractive explanation. The parameters which enter into the neutrino mass matrix 
usually stem from model predictions at high energy scales, such as the scale Mgut 
of grand unification. The measurements and bounds for neutrino masses and lepton 
mixings, on the other hand, determine the parameters at low energy. The high- and 
low-energy parameters are related by the renormalization group (RG) evolution, so that 
low-energy data yield only indirect restrictions for mass models or other high-energy 
mechanisms like leptogenesis [3]. It is well-known that the model independent RG 
evolution between low energy and the lowest see-saw scale can have large effects on the 
leptonic mixing angles and on the mass squared differences, in particular if the neutrinos 
have quasi-degenerate masses |Ul[TllBllUllTn t lll | ll^ U l^ U 14 | ll5 [ ll(i ^ ll7 H 18 ^ ll9 | l2(} H 2 H I22 | l23j. RG 
effects may even serve as an explanation for the discrepancy between the mixings in the 
quark and the lepton sector [21]. 

The RG equations (RGEs) for the neutrino mass operator and for all the other 
parameters of the theory have to be solved simultaneously. The mixing angles, phases 
and mass eigenvalues can then be extracted from the evolved mass matrices. Both 
steps are, however, non-trivial and can only be performed numerically in practice. In 
order to determine the change of the parameters under the RG flow in a qualitative 
and, to a reasonable accuracy, also quantitative way, it is useful to derive analytical 
formulae for the running of the masses, mixing angles and phases. This was done in fTTI] 
assuming CP conservation and in ^T] for the general case. We modify the derivation 
of JT] by a step which simplifies the formulae that arise after explicitly writing out the 
dependence on the mixing parameters. These results are exact, and they make it easier 
to derive simple approximations in the limit of small $13. These approximations are very 
useful in understanding the RG evolution of the phases and the phase dependence of the 
evolution of other parameters. For example, we find that the phases show significant 
running. Consequently, vanishing phases at low energy appear unnatural unless exact 
CP conservation is a boundary condition at high energy, which seems unlikely, since the 
CP phase in the quark sector is sizable. The presence of CP phases at low energies has 
significant impact on observations |2SII2S1I2Z| ■ 

The outline for the paper is: In Sec. El we present analytical formulae for the RG 
evolution of the neutrino masses, leptonic mixing angles and phases, where an expansion 
in the small angle #13 is performed. This leads to very simple and in most cases accurate 
formulae which are compared with numerical results. Sec. HU is devoted to phenomeno- 
logical consequences for leptogenesis, the WMAP bound, the effective neutrino mass 
relevant for neutrinoless double beta decay and precision measurements of #13 and #23- 
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2 RG Evolution of Leptonic Mixing Parameters and 
Neutrino Masses 



In this study, we will focus on neutrino masses which can be described by the lowest- 
dimensional neutrino mass operator compatible with the gauge symmetries of the SM. 
This operator reads in the SM 



A K 9f ^Lc £Cd( t 



e ba <p a + h.c. 



Lb 



(1) 



and in its minimal supersymmetric extension, the MSSM, 



MSSM 



(2) 



K g f has mass dimension —1 and is symmetric under interchange of the generation indices 
/ and g, e is the totally antisymmetric tensor in 2 dimensions, and £^ is the charge conju- 
gate of a lepton doublet, a, b,c,d e {1, 2} are SU(2)l indices. The double-stroke letters 
I and h denote lepton doublets and the up-type Higgs superfield in the MSSM. After 
electroweak (EW) symmetry breaking, a Majorana neutrino mass matrix proportional 
to k emerges as illustrated in Fig. [T] 



Lc EW symmetry 
> 

breaking 



t i 




Figure 1: Vertex from the dimension 5 operator which yields a Majorana mass matrix for the light 
neutrinos. 



The above mass operator provides a rather model-independent way to introduce 
neutrino masses as there are many possibilities to realize it radiatively or at tree-level 
within a renormalizable theory (see e.g. [2H]). The tree-level realizations from integrating 
out heavy singlet fermions and/or Higgs triplets naturally appear for instance in left- 
right-symmetric extensions of the SM or MSSM and are usually referred to as type I 
and type II see-saw mechanisms. 

The energy dependence of the effective neutrino mass matrix below the scale where 
the operator is generated (which we will call Mi in the following) is described by its 
RGE. At the one-loop level, this equation is given by [2313011211122! 

16tt 2 ^ = C {Y}Y e ) T k + Ck (Y}Y e ) + an, (3) 



2 



where t = ln(/j,/ /j, ) and \i is the renormalization scale 1 and where 
C = 1 in the MSSM , 

C = -- in the SM . (4) 

In the SM and in the MSSM, a reads 

«sm = -3g 2 + 2{y 2 + yl + y 2 e ) + 6(y 2 + y 2 + y 2 c +y 2 + y 2 + y 2 u )+X, (5a) 
6 

«mssm = --^g 2 - Qgi + 6 (y 2 + y 2 + vl) • (5b) 

Here Y) (/ 6 {e, d, u}) represent the Yukawa coupling matrices of the charged leptons, 
down- and up-type quarks, respectively, denote the gauge couplings 2 and A the Higgs 
self-coupling in the SM. We work in the basis where Y e is diagonal. 

The parameters of interest are the masses, which are proportional to the eigenvalues 
of k and defined to be non-negative, as well as the mixing angles and physical phases of 
the MNS matrix [HU 

Umns = y(ei 2 ,^3^23,5)diag(e-^/ 2 ,e-^/ 2 ,l), (6) 

which diagonalizes k in this basis. V is the leptonic analogon to the CKM matrix 
in the quark sector. The parametrization we use will be explained in more detail in 
App. El Currently, we learn from experiments that there occur two oscillations with 
mass squared differences Am 2 ol and Am 2 tm and corresponding mixing angles 6*12 and 
823 , respectively. For the third mixing angle #13 and the absolute scale of light neutrino 
masses, there are only upper bounds at the moment (see Tab. [T] for the present status). 



2.1 The Analytical Formulae 

In this section, we present explicit RGEs for the physical parameters. They determine 
the slope of the RG evolution at a given energy scale and thus yield an insight into 
the RG behavior. The derivation will be discussed in App. El Note that a naive linear 
interpolation, i.e. assuming the right-hand sides of the equations to be constant, will 
not always give the correct RG evolution. As we will show later, this is mainly due to 
large changes of 612 and the mass squared differences. In the following, we will neglect 
y e and against y T and introduce the abbreviation 

^ m atm 



1 In the MSSM, the RGE is known at two-loop In this study, we will, however, focus on the 
one-loop equation. 

2 We are using GUT charge normalization for g±. 
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Table 1: Experimental data for the neutrino mixing angles and mass squared differences. For the 
solar angle 812 and the solar mass squared difference, the LMA solution as confirmed by KamLAND 
is shown. The results stem from the analysis [35] of the recent KamLAND and the SNO data, the 
Super-Kamiokande atmospheric data and the CHOOZ experiment pT7] . 



whose LMA best-fit value is about 0.03. In order to keep the expressions short, we 
will only show the leading terms in an expansion in the small angle #13 for the mixing 
parameters. In almost all cases they are sufficient for understanding the features of the 
RG evolution. 3 In all cases except for the running of the Dirac phase 5, the limit #13 — > 
causes no difficulties, the subtleties arising for 5 will be discussed in Sec. 12.4.11 We 
furthermore define rrii(t) := v 2 Ki(t)/A with v = 246 GeV in the SM or v = 246 GeV-sin (3 
in the MSSM and, as usual, Am 2 ol := m\ — m\ and Am 2 tm := mj — m\. Note that 
our formulae cannot be applied if one of the mass squared differences vanishes. For a 
discussion of RG effects in this case, see e.g. [7JIH1III1I221IIIE] • With these conventions, we 
obtain the following analytical expressions for the mixing angles: 

e a = -M^^ K^^I- .^), (8 ) 



'13 



sin 2#i2 sin 29- 



m 3 



x 



Cy 2 

x [mi cos(y?i - 5) — (1 + £) m 2 cos(<^ 2 - 5) - (m 3 cos 5} + 6{Q\?) , (9) 



/03 



Cy 2 
32tt 2 

+ #{013 



sin 29 



28 



c 12 |m 2 e 



1-P2 



+ m 3 \ 2 + s{ 2 



mi e Vl + 777.3 1 

r+c 



(10) 



Note that in order to apply Eq. (0) to the case 6*13 = 0, where S is undefined, the 
analytic continuation of the latter, which will be given in Eq. (|25jh has to be inserted. 

3 The exact formulae, from which we have derived the analytical approximations presented here, can 
be obtained from the web page http://www.ph.tum.de/~mratz/AnalyticFormulae/. 
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The ^(#13) terms in the above RGEs can become important if #13 is not too small and 
in particular if cancellations appear in the leading terms. For example, this is the case 



for \ipi — tp 2 \ = 7T in ([ 
phase is given by 

Cy 2 5^ 



as we will discuss below in more detail. The RGE for the Dirac 
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For the physical Majorana phases, we obtain 
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12a) 
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(13) 



(14) 



We would like to emphasize that the above expressions do not contain expansions in 
(, i.e. their ( dependence is exact. In many cases, they can be further simplified by 
neglecting ( against 1 without losing much accuracy. Note that singularities can appear 
in the £?(0 13 )-terms at points in parameter space where the phases are not well-defined. 
For the masses, the results for y e = y^ = but arbitrary 13 are 



[a + Cy 2 T {2sl 2 s\ 3 + F x )] mi , 
[a + Cy 2 (2c 2 i 2 s 2 23 + F 2 )]m 2 , 



167T nil 
16% 2 fn 2 

16tt 2 m 3 = [a + 2Cy 2 T c 2 13 c 2 23 ] m 3 



(15a) 
(15b) 
(15c) 
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where F\ and F 2 contain terms proportional to sin #13, 

Fi = —S13 sin 26*12 sin 26*23 cos 5 + 2s^ 3 c^ 2 C2 3 , (16a) 
F 2 = S13 sin 26*12 sin 26*23 cos 5 + 2s^ 3 C23 • (16b) 

These formulae can be translated into RGEs for the mass squared differences, 

8vr 2 ^Am2 ol = a Am 2 sol + Cjfi [2s 2 23 (mj c\ 2 - m\ s\ 2 ) + F sol ] , (17a) 

87(2 Jt Am ' tm = Am ' tm + ^ ^ " ^2 C 12 S 23 + ^atm] , (17b) 

where 

Fsoi — { m i + m 2 ) s n sin 2^12 sin 26*23 cos 5 

+ 2s 13 C 23 ( m 2 S 12 - m l C 12) , ( 18a ) 

-^atm = — m 2 S13 sin26* 12 sin 26*23 cos 5 — 2m\ s\ 3 s\ 2 (? 23 . (18b) 

2.2 Generic Enhancement and Suppression Factors 

From Eqs. (|8*|)— (|14|) it follows that there are generic enhancement and suppression factors 
for the RG evolution of the mixing parameters, depending on whether the mass scheme 
is hierarchical, partially degenerate or nearly degenerate. We have listed these factors 
in the approximation of small 6*13 in Tab. El They can be compensated by cancellations 
due to a special alignment of the phases. For example, an opposite CP parity of the 
first and second mass eigenstate, i.e. \<pi — (p 2 \ = 7T, results in a maximal suppression 
of the running of the solar mixing angle, which has been pointed out earlier in papers 
like [TT } I39 | I13 | I17|. Nevertheless, Tab. |2] allows to determine which angles or phases have 
a potential for a strong RG evolution. Obviously, the expressions for S are not applicable 
for #13 = 0. This special case will be discussed at the end of Sec. 12.4.11 

Let us consider some numerical values in order to estimate the size of RG effects. 
The SM r Yukawa coupling is yf M = ^m T fa 0.01. Thus, the typical factor in the 
formulae for the mixing angles and phases amounts to 

Ut -0.5-10- 6 . (19) 



64tt 2 

In the MSSM it changes to 

-^^0.3-KT 6 (l+tan 2 /3) . (20) 

If the running was purely logarithmic, it would yield a factor of 
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Tabic 2: Generic enhancement and suppression factors for the RG evolution of the mixing parameters. 
A '1' indicates that there is no generic enhancement or suppression, 'n.h.' and 'p.d.(n.)' denote the 
hierarchical and partially degenerate mass spectrum in the case of a normal hierarchy, i.e. ml <C Am^ ol 
or Ar7ig ol <C m\ < Am^tm- 'i-h.' and 'p.d.(L)' denote the analogous spectra in the inverted case, i.e. 
m§ <€; Am.g ol or Arn^, <§; m\ < Am^ tm . Finally, 'd.' means nearly degenerate masses, Am1 tm <C ml ~ 

2 2 2 

777,2 ~ m 3 ~ TO • 

for Mi = 10 13 GeV. If we assume that the solar and atmospheric angle are large and 
that the phases do not cause excessive cancellations, then multiplying the above two 
contributions with the enhancement factor r en h from Tab. |2] yields a rough estimate for 
the change of the angles and phases due to the RG evolution, 

A RG ~10- 5 (l+tan 2 /3)r enh . (22) 

Of course the factor 1 + tan 2 (3 has to be omitted in the SM. It is immediately clear that 
even in the MSSM with very large tan /3 no significant change occurs if the enhancement 
factor is 1 or less - except maybe for 9 X3 , where even a change by 1° could be interesting. 
However, for quasi-degenerate neutrinos large enhancement factors are possible. As an 
example, let us estimate the size of the absolute neutrino mass scale (the 'amount of 
degeneracy') needed for a sizable RG change of #12, say 0.1 ~ 6°. In the SM, this requires 
T e nh ~ 10 4 , corresponding to a neutrino mass of the order of leV, which is excluded 
by WMAP and double beta decay experiments. On the other hand, in the MSSM this 
mass scale can easily be lowered to about 0.1 eV with tan/3 as small as 8. 

2.3 Discussion and Comparison with Numerical Results 

We now study in detail the running of the mixing angles and masses, in particular the 
influence of the phases. The RG evolution of the phases will be studied separately in 
Sec. 12.41 We solve the RGEs for the neutrino mass operator and for the other parameters 
numerically and compare the results with those obtained from the analytical formulae 
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of Sec. 12.11 For the numerics we follow the 'run and diagonalize' procedure, i.e. we first 
compute the running of the mass matrix and then extract the evolving mass eigenvalues 
and mixing parameters. The algorithm used for this is described in App. [X] As an 
example, we consider the MSSM with tan/3 = 50, a normal mass hierarchy for the 
neutrinos, mi = 0.1 eV for the mass of the lightest neutrino, and a mass of about 
120 GeV for the light Higgs. These boundary conditions are given at the electroweak 
scale, i.e. we calculate the evolution from low to high energies. Below the SUSY-breaking 
scale, which we take to be 1.5 TeV, we assume the SM to be valid as an effective theory 
and use the corresponding RGEs. Above, we apply the ones of the MSSM. 



2.3.1 RG Evolution of 12 

From Tab. |21 we see that the solar angle 0i2 generically has the strongest RG effects 
among the mixing angles. The reason for this is the smallness of the solar mass squared 
difference associated with it, in particular compared to the atmospheric one, which 
leads to an enhanced running for quasi-degenerate neutrinos and for the case of an 
inverted mass hierarchy. Furthermore, it is known that in the MSSM the solar angle 
always increases when running down from Mi for 13 = 20 . This is confirmed by 
our formula (jSJ). From the term \mie lipi + rri2e lip2 \ 2 in Eq. (JBJ), we see that a non-zero 
value of the difference \(pi — y^l of the Majorana phases damps the RG evolution. The 
damping becomes maximal if this difference equals tt, which corresponds to an opposite 
CP parity of the mass eigenstates mi and m 2 . This is in agreement with earlier studies, 

e.g. imESUEllIZI. 

Let us now compare the analytical approximation for 812 of Eq. (JHJ) with the numerical 
solution for the running in the case of nearly degenerate masses, which is shown in Fig. El 
in detail. The dark-gray region shows the evolution with LMA best-fit values for the 
neutrino parameters, #13 varying in the interval [0°,9°] and all CP phases equal to 
zero. The medium-gray regions show the evolution for \ipi — ip 2 \ G {0°, 90°, 180°, 270°}, 
013 e [0°,9°] and 5 G {0°, 90°, 180°, 270°}, confirming the expectation of the damping 
influence of (pi and <p2- The flat line at low energy stems from the SM running below 
^susy, which is negligible as we have seen earlier. Note that the numerics never yield 
negative values of 612 due to the algorithm used for extracting the mixing parameters 
from the MNS matrix, which guarantees < 6i 2 < 45° (see App. lA.3l for further details). 

As can be seen from the relatively broad dark-gray band in the figure, the ^(#13)- 
term in the RGE is quite important here. The dominant part of this term is 

Cy 2 m 2 + mi . <£i - <p 2 
T = sm 2023 cos = x 

6ZTi z TO2 — Tn% 2 

x (^cos2#i 2 coso cos hsmo sin J • #13 . [26) 

Clearly, the RG evolution of 6*12 is independent of the Dirac phase 5 only in the ap- 
proximation 0i3 = 0. The largest running, where 0i2 can even become zero, occurs for 
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Figure 2: RG evolution of Q\i in the MSSM with tan/3 = 50, a normal mass hierarchy and m\ = 
0.1 eV. The dark-gray region shows the evolution with best-fit values for the neutrino parameters, 
#13 € [0°,9°] and all CP phases equal to zero. The medium-gray regions show the evolution for 
\(pi - ip 2 \ = 0°, \cpi - tp 2 \ £ {90°, 270°} and |</?i - tp 2 \ = 180°. They emerge from varying 13 e [0°,9°] 
and 6 £ {0°, 90°, 180°, 270°}. The light-gray regions can be reached by choosing specific values for 
the CP phases different from the ones listed above. The dashed line shows the RG evolution with 
\<Pi — ¥>2 1 = 0, 6>i3 = 9° and S = 180°. Note that for the numerics we use the convention where 812 is 
restricted to the interval [0°,45°], so that the angle increases again after reaching 0. The dotted line 
shows the evolution with \<pi — (p 2 \ = 90° and #13 = 0°. 

9 13 as large as possible (9°), 5 = tc and <p% — ip 2 = 0. In this case the leading and the 
next-to-leading term add up constructively. It is also interesting to observe that due to 
^(#13) effects 9 12 can run to slightly larger values. The damping due to the Majorana 
phases is maximal in this case, which almost eliminates the leading term. Then, all the 
running comes from the next-to- leading term (|23j). 

In the inverted scheme, mi 3> m 2 — mi always holds, so that large RG effects are 
generic, i.e. always present except for the case of cancellations due to Majorana phases. 
For a normal mass hierarchy with a small mi, the running of the solar mixing is of 
course rather insignificant. 

Finally, we would like to emphasize that it is not appropriate to assume the right- 
hand sides of Eq. (jHJ) and Eq. (f23j) to be constant in order to interpolate 9 12 up to a high 
energy scale, since non-linear effects especially from the running of sin 2#i2 and Am^ 
cannot be neglected here. This is easily seen from the curved lines in Fig. 121 

2.3.2 RG Evolution of 13 

The analytical approximation for #13 is given in Eq. (jHJ). As already pointed out, in order 
to apply it to the case #13 = 0, where 5 is undefined, the analytic continuation of the 
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latter has to be inserted. It will be given in Eq. (|2l)j) in section 12. 4. 11 where the phases 
are treated in detail. The comparison with the numerical results in Fig. |3] shows that 
above Msusy the angle runs linearly on a logarithmic scale to a good approximation. 
Thus, using Eq. © with a constant right-hand side yields pretty accurate results. With 
if i 7^ ip2, significant RG effects can be expected for nearly degenerate masses. This is 
confirmed by the light-gray region in Fig. El 

The fastest running occurs if ipi — <p% = n and ipi — 8 G {0,7r}, so that the terms 
proportional to mi and in the RGE are maximal and add up. Interestingly, cancel- 
lations between the first two terms in the second line of Eq. (j^J) appear for (pi = <p2, 
in particular if all phases are zero. If so, the leading contribution to the evolution of 
#13 is suppressed by an additional factor of £. This suppression is in agreement with 
earlier studies, for instance [123131], where it was discussed for the CP-conserving case 
(pi = (f2 = 7r, which implies an opposite CP parity of m$ compared to the other two 
mass eigenvalues. Such cancellations cannot occur for a strong normal mass hierarchy, 
since then the evolution is dominated by the term proportional to m 2 in Eq. (JHj). 

Besides, #13 runs towards smaller values in the MSSM with zero phases and a normal 
hierarchy, because mi < m 2 , so that the second line of the RGE is negative. This yields 
the dark-gray region in Fig. El 4 As #13 can always be made positive by a suitable 
redefinition of parameters, the sign of #13 is irrelevant for 813 = 0. 

For an inverted hierarchy, the situation is reversed, since Am 2 tm is negative then. 
For a small m.3, the running is highly suppressed in this case, because the leading term 
is proportional to m 3 . Then the dominant contribution comes from the ^(fl^-term 
unless 8 13 is very small as well. 

Future experiments will probably be able to probe sin 2 2#i3 down to 10~ 4 , corre- 
sponding to #i3 ~ 5 • 1CT 3 ~ 0.3°. Consequently, even RG changes of this order of 
magnitude could be important, since a low-energy value smaller than the RG change 
would appear unnatural. This will be discussed in more detail in Sec. 13.31 

2.3.3 RG Evolution of 23 

The analytical RGE for 823 can be found in Eq. ([10)1 . Again, the comparison with the 
numerical results (see Fig. OJ) shows that to a good approximation the angle runs linearly 
on a logarithmic scale above Msusy- The sign of Am 2 tm is very important here. For a 
normal mass spectrum, the leading term is always negative in the MSSM, so that 823 
decreases with increasing energy, while for an inverse spectrum the situation is exactly 
reversed, so that 823 becomes larger than 45° if one starts with the LMA best-fit value 
at low energy. 

From Eq. (|TU|) we expect that switching on the phases ipi and <f2 always reduces the 
running of 823 for nearly degenerate masses. This is confirmed by the light-gray region 

4 The relatively large slope of its upper boundary is due to the ^(6*13) contribution to the RGE. 
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Figure 3: RG evolution of #13 and 6*23 in the MSSM with tan/3 = 50, a normal mass hierarchy and 
mi =0.1 eV. The dark-gray region shows the evolution with best-fit values for the neutrino parameters, 
013 € [0°,9°] and all CP phases equal to zero. For the #23 case, we just obtain a thick gray line at the 
bottom of the gray region. The light-gray regions show the evolution, which is possible, if arbitrary CP 
phases are allowed. 

in Fig. El The damping is much less severe for a hierarchical mass spectrum, since either 
mi and m2 or 777.3 are very small then. However, in these cases the running is generally 
expected to be rather insignificant, since according to Tab. Elthe enhancement factor is 
only 1. 

2.3.4 RG Evolution of the Neutrino Mass Eigenvalues 

The running of the mass eigenvalues is significant even in the SM or for strongly hier- 
archical neutrino masses due to the factor a in the RGEs (|15j). Clearly, the evolution 
is not directly dependent on the Majorana phases [H]. This can be understood from 
Eqs. f|B.13jl and (|B.19J) . which show that only the moduli of the elements of the MNS 
matrix enter into rhj. Besides, rh.3 does not depend on 5, since only the moduli of the 
elements of the third column of the MNS matrix are relevant in this case. Of course, 
there is an indirect dependence on the phases, as these influence the running of the 
mixing angles. 

Apart from the MSSM with large tan/5, the running of the mass eigenvalues is 
virtually independent of the mixing parameters, since a is usually much larger than y\. 
In the SM, the Higgs mass influences the running via the self-coupling A - the heavier 
the Higgs, the larger the RG effects. Thus, except for large tan/3 in the MSSM, the 
running is given by a common scaling of the mass eigenvalues , which is obtained by 
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neglecting y T and integrating Eq. (fI5j). 



mj(i) 



exp 



16tt 5 



dr air) 



to 



mi(t ) =: s(t,t )mi(t ) 



(24) 



We plot s in the SM and in the MSSM for various parameter combinations in Fig. 0] 
The three SM curves correspond to different Higgs masses in the current experimentally 
allowed region at 95% confidence level, 114 GeV < m H < 200 GeV |H|. m H = 180 GeV 
is the value for which the self-coupling A stays perturbative up to 10 16 GeV, i.e. A < 1, 
and rriH = 165 GeV is the minimal mass for which A is positive up to 10 16 GeV, so 
that the vacuum is stable in this region (see e.g. [HJI32])- 5 ^ e MSSM, we choose 
ma — 120 GeV for the light Higgs mass, since the allowed range is further restricted 
by the upper limit at about 130 GeV here, and since it influences the evolution of 
the RG scaling only marginally as long as Msusy and Mz differ only by a few orders 
of magnitude. Moreover, further uncertainties due to threshold corrections and the 
unknown value of the SUSY-breaking scale can be equally important as the one due to 
the unknown Higgs mass. The RG enhancement of the masses is smallest if tan/3 pa 10. 
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Figure 4: Scaling of the masses under the renormalization group in the SM and MSSM. The mixing 
parameters are chosen to be the LMA best- fit values (cf. Tab.^), but they influence the running only 
marginally. We further used a SUSY-breaking scale Msusy = 1 TeV. The upper curves show the 
evolution in the SM for run = 114 GeV, run = 165 GeV and ran — 180 GeV, the lower ones correspond 
to the MSSM for tan/3 = 10 and tan/3 = 50 with ran — 120 GeV. These plots apply for all mass 
eigenvalues, except for large tan/3 in the MSSM where the scaling of m.3 is shown (using zero phases). 
Note also that a different SUSY-breaking scale changes the scaling factor in the MSSM. 



5 In some models (see, e.g. 03] for a viable model) A can be larger, in particular if M\ <C 10 16 GeV. 
A negative value of A at high energy implies a metastable vacuum. 
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As already mentioned, substantial deviations from the common scaling arise in the 
MSSM for large tan (3. There is a plethora of effects which can be understood with the 
aid of (fT5|) and ()17|) . In order to give an interesting example, we show the evolution of 
the mass eigenvalues for m m j n = 0.19 eV (where m m j n = min{mi, m2, m^]) in the MSSM 
with tan (3 = 50 in Fig. A particular interesting effect is that for an inverted mass 
spectrum the property |Am^ tm | > Am^ possibly does not survive the RG evolution. 
In other words, what looks like a normal mass hierarchy at high energies turns out to 
become an inverted hierarchy at low energies (cf. Fig. |5(b)| . From the dependence on 
the yl terms (cf. Eqs. (fi7)|) and (fT%|) ). we find that this effect can disappear if S is large. 




2 4 6 8 10 12 2 4 6 8 10 12 

log 10 (///GeV) log 10 (///GeV) 

(a) Normal mass hierarchy (b) Inverted mass hierarchy 



Figure 5: Running of the light neutrino masses for a normal and an inverted mass hierarchy and 
Ti m i n = 0.19 eV in the MSSM with tan/3 = 50 and Mgusy = 1 TeV. The mixing parameters are 
chosen to be the LMA best-fit values. The phases are zero in this example. In the inverted case, Am s 2 d 
becomes greater than |Am| tm |. 

2.3.5 RG Evolution of Am^, 

The RGE for the solar mass squared difference is given in Eq. ()17a|) . In the SM and the 
MSSM with small tan/?, the running is due to the common scaling of the masses de- 
scribed in the previous section and thus virtually independent of the mixing parameters. 
For large tan (3 and nearly degenerate masses, the influence of CP phases, in particular 
the Dirac phase, is crucial. The numerical example in Fig. |U] confirms this expectation 
and furthermore shows that Am;: ol runs dramatically. On the one hand, it can grow by 
more than an order of magnitude. As we have seen in Fig. A?rig 0l can even get larger 
than |Am^ tm |. On the other hand, it can run to at energy scales slightly beyond the 
maximum of 10 13 GeV shown in the figure. For large tan (3, Am^ -C ml and not too 
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small #13, the first term in F so i is essential for understanding these effects, since it is 
proportional to the sum of the masses squared rather than the difference. For 5 = tc and 
#13 near the CHOOZ bound, its sign is negative and its absolute value maximal, which 
causes the evolution of Am^ towards zero. For 5 = 0, the sign becomes positive, so 
that the running towards larger values is enhanced, which explains the upper boundary 
of the light-gray region in Fig. |U1 




2 4 6 8 10 12 

log 10 Gu/GeV) 



Figure 6: RG evolution of Am^ ol in the MSSM with tan/3 = 50, a normal mass hierarchy and 
mi = 0.1 cV. The dark-gray region shows the evolution with LMA best-fit values for the neutrino 
parameters, #13 S [0°,9°] and all CP phases equal to zero. The light-gray regions show the evolution, 
which is possible, if arbitrary CP phases are allowed. 



2.3.6 RG Evolution of Am^ tm 

From the numerical example in Fig. [71 we see that Am^ tm can be damped by the phases, 
but not significantly enhanced. Depending on the CP phases, Am^ tm grows by about 
50% - 95%. Analogously to above, the maximal damping is mainly due to the first term 
in F atm , so that it occurs for large #13 and 5 = 0. Compared to the case of the solar mass 
squared difference, the influence of 5 is generically smaller here, because Aml tm /mj is 
larger and because the phase-independent terms in the RGE do not nearly cancel. 
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Figure 7: RG evolution of Am^ tm in the MSSM with the same input parameters as in Fig. 



2.4 RG Running of the Dirac and Majorana Phases 

Most earlier studies of RG effects either neglected phases or concentrated on the special 
case of a Majorana parity, where one or both of the Majorana phases are n. We have 
seen that they can have a dramatic influence on the running of the masses and mixings. 
Moreover, many effects are affected by phases, e.g. neutrinoless double beta decay, or 
require phases, e.g. leptogenesis. 6 

Of course, if the phases are given at some scale, they also change due to the RG 
evolution. We now discuss the running of the phases themselves and give numerical 
examples. In general, a significant evolution of the phases is expected for nearly de- 
generate and inverted hierarchical mass patterns, since the RGEs ()lip -([13| ) contain the 
ratios m^j/Am^. 

2.4.1 RG Evolution of the Dirac Phase 

The running of the Dirac phase 5 is given by Eq. (jllj) for y e = = 0. An interesting 
possibility is the radiative generation of a Dirac phase by Majorana phases [H]: A 
non-zero 6 is produced by RG effects, since some of the terms in the RGE (fTTj) do not 
vanish for 8 — > 0. Fig. |H] shows an example. The most important term in this context 

6 Clearly, the phases relevant for leptogenesis are those of the 'right-handed' sector and therefore in 
general not directly related to the phases considered here |44II45| . However, as the left-handed sector 
with its - in principle - observable phases is related to the right-handed one by the see-saw relation, it 
is reasonable to assume that non- vanishing right-handed phases imply non-zero <5, tpi and/or if2- An 
explicit relation which supports this point of view is specified in, e.g., |46| . 
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is the first one in 5^°\ As it is proportional to sin(</? 1 — </? 2 ), the effect is suppressed for 
ifi = ip 2 . For small but non-zero values of #13 , the term involving S^~^ also contributes 
significantly because of the factor 9^. For ipi = ip 2 , this contribution is suppressed as 
well, since the parts proportional to mi and m 2 , respectively, nearly cancel. 
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Figure 8: Radiative generation of a Dirac phase in the MSSM with tan/3 = 30 and a normal hierarchy. 
Here the running is from high to low energy, i.e. the boundary conditions are given at the see-saw scale. 
S is zero there but large at Mz- The other starting values are 9i 2 = 18°, 6*i 3 e {1°,3°,6 }, #23 = 34°, 
mi = 0.17 eV, Am^ m = 3.8 • 10- 3 eV 2 , Ara 2 o] = 5.7- 10~ 4 cV 2 , ipx = 16°, ip 2 = 140°. 

In the case of an inverted hierarchy with tan/3 varying between 30 and 50, Dirac 
phases of about 15° to 30° can be generated. Now the term involving S^ 1 ^ receives 
an additional suppression from the small value of m 3 , so that the subleading effects 
described above become unimportant. Hence, the running of 5 is independent of #13 and 
depends only on the difference of the Majorana phases to a very good approximation. 

Before we turn to the evolution of the Majorana phases, let us discuss some further 
properties of the RGE for 5 that are also valid beyond the special case of a radiative 
generation of this phase. To start with, the most important term in 5 depends only 
on the difference of the Majorana phases. Consequently, the evolution is expected to 
stay roughly the same if both phases change by the same value. A comparison with 
numerical results shows that this is true only to a first approximation. If one starts with 
</?2 = and increments it step by step, the running of 5 is increasingly damped. The 
main reason for this is the second term in square brackets in (the one proportional 
to m 2 ), whose sign is opposite to that of the leading term for 5 < ip 2 . This term grows 
with </? 2 , while the previous one (proportional to mi) does not change much as long as 
ipi is close to 90°. The situation can be very different for smaller values of #13. Now the 
initial rise of 5 is enhanced, so that it can become larger than (p 2 . Then the sign of the 
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aforementioned second term in square brackets changes, so that it no longer damps the 
evolution but amplifies it. 

With a strong normal hierarchy, RG effects are usually tiny. The running of the 
Dirac phase is one of the few examples where this is not always the case. Due to the 
terms proportional to d^ 3 in the RGE, a significant evolution is possible for small #13. 
However, one has to keep in mind that a measurement of 5 is very hard in this case. 

Regardless of the mass hierarchy, the limit 6*13 — > is dangerous, because in this 
case the RGE fjlljl diverges. However, we can show that S remains well-defined: The 
derivative of the MNS matrix U is given by ()B.9j) . U = U ■ T, where U and T are 
continuous. Hence, Uis(t) describes a continuously differentiate curve in the complex 
plane. Consequently, # 13 and 6 are continuously differentiate even for # 13 = 0, if 5 
is extended continuously at this point. Note that restricting the parameters to certain 
ranges can nevertheless result in discontinuities. For example, if the RG evolution causes 
#13 to change its sign and if we demand < #13 < |, then there will be a kink in the 
evolution of #13 and 5 will jump by it. However, even in the presence of such artificial 
discontinuities there must still be finite one-sided limits for 5 and 5 as # 13 approaches 0. 

The limit for 5 is determined by the requirement that 8 remains finite. Then the 
divergence of 9^ has to be canceled by f^ -1 ). For ipi = ip 2 = 0, this obviously implies 
5 = or 5 = 7r. In the general short calculation yields 

m 1 cos(p 1 - (l + ()m 2 cosip 2 -(m 3 

cotd = : — ; . (25) 

rrii sin (pi — (1 + (J ?»2 sin (p 2 

Due to the periodicity of cot, there are two solutions differing by tt, corresponding to 
the different limits on the two sides of a node of 6 13 . 



2.4.2 RG Evolution of the Major ana Phases 

While the RGEs for the Majorana phases are somewhat lengthy, there is a simple ex- 
pression for the running of their difference for small #13, 

0i - 02 = -r~k ~\ — 2~ cos 26l i2 sin2 #23 sin(v9i - v? 2 ) + ^(#13) • (26) 
4vr 2 Am 2 sol 

It shows that for 813 = 0, the phases remain equal, if they are equal at some scale. 
Obviously, (p\ — up 2 > for (pi > (p 2 and vice versa, which means that the difference 
between the phases tends to increase with increasing energy. In other words, a large 
difference at the see-saw scale becomes smaller at low energy. An example is shown in 
Fig.0 

If (fi — if 2 is not too small, a non-zero 613 tends to damp its running. This is due to a 
term in the RGE for ipi whose sign is opposite to that of the leading one in Eq. ()26|) and 
which is proportional to sin 613 cot #12. This term can grow important if 6i 2 becomes 
small with increasing energy. 
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Figure 9: Running of the Majorana phases in the MSSM with a normal hierarchy, tan (3 = 50, ip% = 75°, 
if 2 = 70°, #i3 = 0, mi = 0.15 eV, and LMA best-fit values for the mass squared differences, 812 and 623 
at Mz- RG effects are substantial, and the difference ipi — ip2 increases with increasing energy. 

For (fx = cp2 the evolution of the Majorana phases is suppressed, since the leading 
terms in the RGEs (j!3|) and (fT3jl are zero then. However, for larger tan f3 RG effects are 
still important. Non-linear effects caused by the decrease of the solar and atmospheric 
mixing angles are essential here, as the initial slope of the curves is extremely small due 
to the suppression by sin #13 and cos26 , 2 3. For 9i 3 = 5°, the second line in the RGE 
and the terms proportional to sin# 13 are about equally important for the running of tp\. 
The evolution of (p 2 is virtually independent of 6^3, since the respective terms are not 
multiplied by cot 612, which again can become large as the energy increases because of 
the diminishing 812, but by tan #12, which remains smaller than 1. 

In principle, it is also possible to generate Majorana phases radiatively, if the CP 
phase is non-zero. However, it follows from the discussion in the previous paragraph 
that this only happens via terms proportional to sin^. 

3 Some Applications 

The discussed RG effects obviously have important implications whenever masses and 
mixings at different energy scales enter the analysis. 
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3.1 Relating the Leptogenesis Parameters to Observations 

One of the most attractive mechanisms for explaining the observed baryon asymmetry 
of the universe, rjs = (6.5l ; 8 ) • 1CT 10 |H], is leptogenesis [Sj. In this scenario, t]b is 
generated by the out-of-equilibrium decay of the same heavy singlet neutrinos which are 
responsible for the suppression of light neutrino masses in the see-saw mechanism. The 
masses of the heavy neutrinos are typically assumed to be some orders of magnitude 
below the GUT scale. 

Though the parameters entering the leptogenesis mechanism cannot be completely 
expressed in terms of low-energy neutrino mass parameters, it is possible to derive 
bounds on the neutrino mass scale from the requirement of a successful leptogenesis [48J. 
Since, as we demonstrated in Sec. 12.3.41 the neutrino masses experience corrections of 
about 20-25% in the MSSM or more than 60% in the SM, we expect the corrections for 
such bounds to be sizable. 

The maximal baryon asymmetry generated in the thermal version of this scenario is 
given by [131231111 

r/ max ~ 0.96 ■ 10~ 2 ef 3 * Kf . (27) 

Kf is a dilution factor which can be computed from a set of coupled Boltzmann equations 
(see, e.g. [SI])- In |48j . an analytic expression for the maximal relevant CP asymmetry 
was derived, 



max/ ~ \ 3 M 1 m 3 

£r(mi,m 3 ,mi) 



16tt ( v /y/2) 
which refines the older bound 



2 2\ 1/2' 

1 mi ( I _|_ m 3 ~~ m l 



m 3 \ m\ 



(28) 



<TK«J = ~ % + Am "' (29) 

167T (v/V2) 2 m 3 

and is valid for a normal mass hierarchy in the SM as well as in the MSSM. 7 rh\ is 
defined by 

fi, = (30) 

with wd ~ Y v being the neutrino Dirac mass and typically lies between rri\ and m 3 . It 
can be constrained by the requirement of successful leptogenesis because it controls the 
dilution of the generated asymmetry. The authors of f^HJ introduced the 'neutrino mass 
window for baryogenesis' which corresponds to the region in the rh\-M\ plane allowing 
for successful thermal leptogenesis. The shape and size of the 'mass window' depends 

7 To use these formulae in our conventions for the inverted scheme, one would have to replace 
(mi,m 2l m 3 ) -> (m 3 ,m 1 ,m 2 ). 
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on m = \Jm\ + m\ + mf, i.e. it becomes smaller for increasing m, and fn > 0.2 eV is 
not compatible with thermal leptogenesis. 

The calculations relevant for leptogenesis, however, refer to processes at very high 
energies, and therefore the RG evolution of the input parameters has to be taken into 
account [S2] • The correct procedure would be to assume specific values for the neutrino 
mass parameters at low energy, taking into account the experimental input, evolve them 
to the scale Mi and test the leptogenesis mechanism using these values. As the full 
calculation is beyond the scope of this paper, we present the evolution of the relevant 
mass parameters, i.e. the light neutrino masses, to the leptogenesis scale Mi and estimate 
the size of the error arising if RG effects are neglected. 

As discussed in Sec. 12.3.41 there are basically two cases which have to be distin- 
guished, the case of the SM or the MSSM with small tan (3, and the case of the MSSM 
with large tan/?. 

In the first case, running effects can be understood to arise due to the rescaling of 
the light neutrino mass eigenvalues under the renormalization group. From Eq. (|2l?|) it 
is clear that the maximal CP asymmetry scales like the masses. This statement also 
holds for the asymmetry from Eq. ()28|) . if rh\ is a linear combination of the light mass 
eigenvalues. Hence, the RG yields an enhancement of the CP asymmetry of between 
10% and 80%, which can be read off from Fig. |3J These effects are almost completely 
independent of the low-energy CP phases. On the other hand, the dilution factor Kf is 
expected to become tiny since larger mass eigenvalues imply larger Yukawa couplings, 
which makes the washout more efficient. This expectation is substantiated by the fact 
that fn, which controls an important class of washout processes, also increases under the 
renormalization group, i.e. it scales like the masses. As a detailed numerical calculation 
of the dilution factor is beyond the scope of this paper, we refer to jHI], from which we see 
that in the region of interest, i.e. the edge of the mass window, Kf decreases exponentially. 
From this behavior, which is also in accordance with the analytic approximations (see, 
e.g. [SSlElj), we expect that the neutrino mass window for baryogenesis will rather 
shrink than become larger when RG effects are properly taken into account. 

In the second case, i.e. in the MSSM for large tan /3, we distinguish between hierar- 
chical and degenerate mass spectra. In the hierarchical spectrum, the running of e™ 8 * 
is to a high accuracy given by the running of m^, 8 so that in this case Fig. 0] yields the 
relevant plot. The scaling depends on tan/3. In order to illustrate this dependence, we 
pick M 1 = 10 10 GeV and plot fn Tcl := m(10 10 GeV)/m(M z ) in Fig. |10(a)| as a function 
of tan/3, including small values of this parameter as well. It is clear that fn ~ so 
that Fig. |10(a")1 also shows the scaling of e™ ax . Since tan/5 = 10 and tan/5 = 50 corre- 
spond to extreme cases, the scaling factor for different Mi can be read off from Fig. 0] 
by interpolation. 

8 For an inverted hierarchy, mi has to be used instead, whose evolution is approximately the same 
as that of to 3 here. 
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In the case of a quasi-degenerate mass spectrum (and large tan/5), the CP asym- 
metry can run stronger than the average mass scale because, as we already have seen 
in Sec. 12.3.51 and I2.3.6[ the mass squared differences can experience a stronger RG 
enhancement than the squares of the mass eigenvalues. We show the evolution of 
e iel ■= ef^lO 10 GeV)/ef ax (M z ) in Fig. [lO(b)| To produce this plot, we employed 
()29p and inserted the running mass parameters. For this combination of parameters, 
the low-energy phases do influence the evolution of e re \ by damping its running, and the 
plot shows the maximal evolution, which means that the phases are simply set to zero. 
The running effects are even larger for the new bound ()28)1 . since it is more sensitive 
to the mass splittings than the old one. More precisely, for highly degenerate mass 
spectra it is much smaller than the old one and the degeneracy can be lifted by running 
effects. This strong enhancement of the CP asymmetry may even overcompensate the 
decrease of the dilution factor for large tan (3, so that the parameter region compatible 
with thermal leptogenesis grows. 
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(a) tan /3-dependence of the scaling of (b) mi-dependence of e IC \ for 

to. tan/3 = 50. 



Figure 10: Radiative enhancement of to and the CP asymmetry in the MSSM. We show only the case of 
a normal mass hierarchy, since an inverted hierarchy yields virtually the same plot. We assume a SUSY- 
breaking scale Msusy = 1 TeV, a leptogenesis scale of 10 10 GeV, and zero phases. The mixing angles 
and mass squared differences are the LMA best-fit values. We define e rc i := e™ ax (10 10 GeV) / ' e™ a *(Mz)- 
In the case of degenerate masses (see the right part of plot (b)), e™ ax can run stronger than the mass 
eigenvalues since the mass squared differences can have a stronger dependence on the renormalization 
scale than the squares of the mass eigenvalues (cf. Fig. E}. 



Altogether, we have presented the relevant mass parameters at the scale of lepto- 
genesis, thus making it convenient to take into account RG effects in future studies. 
Moreover, we have estimated the impact of the renormalization effects, and found that 
there are two effects in opposite directions: The CP asymmetry is enhanced because the 
mass squared differences increase, and the dilution of the baryon asymmetry is more 



21 



effective since the overall mass scale rises due to RG effects. As the dependence of the 
dilution factor on the mass scale is stronger than that of the CP asymmetry, we expect 
the mass window for baryogenesis to shrink when RG effects are included in the analysis. 
An exception is the case of large tan (3, where the situation is more complicated. 

Note also that there exist different, non-thermal baryogenesis mechanisms jHS] in 
which the masses of the light neutrinos may be almost degenerate [56]. In these kinds 
of scenarios, RG effects increase the baryon asymmetry, since £\ increases, while the 
effects from the expected decrease of the dilution factor do not occur. 

3.2 RG Evolution of Bounds on the Neutrino Mass Scale 

The absolute neutrino mass scale at low energy is restricted by low-energy experiments 
such as searches for Qvfifi decay and cosmo logical observations. As usual, the RG evo- 
lution of the results has to be taken into account in order to translate the experimental 
results into constraints on high-energy theories. 

3.2.1 Neutrinoless Double Beta Decay 

The amplitude of 0^/3/9 decay is proportional to the effective neutrino mass 
(m„) = (m„)n = l^t/^m; 

i 

= \m x c\ 2 c\ 3 e 1 ^ +m 2 s 2 12 cl 3 e i(p2 + m 3 s 2 13 e 2iS | , (31) 

where U is the MNS matrix. Instead of inserting the lengthy RGEs for all the quantities 
in the second line in order to calculate the RG evolution of (m v ), it is much more 
convenient to use Eq. Q, which directly yields 

16tt 2 A (m v ) = {2Cy 2 e + a) (m„> . (32) 

As the first term is negligible, the RG change of the effective neutrino mass is basically 
caused by the universal rescaling of the neutrino masses alone. It is completely inde- 
pendent of the other neutrino mass parameters, since neither the running of y e nor that 
of the terms in a is sensitive to them. Besides, the value of tan j3 is not very important 
here, because y 2 is always tiny and a contains only the up-type quark Yukawa couplings 
in the MSSM. However, there is a dependence on the Higgs mass in the SM. 

Currently, the best experimental upper limit on the effective neutrino mass is about 
(m v ) < 0.35 eV jSZlEHI, with some uncertainty due to nuclear matrix elements. Fig. HP 
shows the running of this limit in the SM and the MSSM. As it is very close to the 
best-fit value of the recently claimed evidence for double beta decay, (m v ) = 0.39 eV 
the evolution of the latter is nearly identical. The SM plot contains three curves 
corresponding to different Higgs masses in the current experimentally allowed region. In 
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the MSSM, the light Higgs mass is chosen to be about 120 GeV. The running is much 
more significant in the SM than in the MSSM because of the contribution of the Higgs 
self-coupling. 
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(a) SM (b) MSSM (tan/3 = 50, 

Msusy = 1.5 TeV) 

Figure 11: Extrapolation of the experimental upper limit on the effective neutrino mass for Qvf3(3 decay, 
(m u ) — 0.35 eV, to higher energies. The SM curves correspond to Higgs masses of 114 GeV, 165 GeV 
and 190 GeV (from bottom to top). In the MSSM, a light Higgs mass of 120 GeV is used. 



3.2.2 WMAP Bound 

Combining the observations of the cosmic microwave background by the WMAP satellite 
with other astronomical data allows to place an upper bound of about 0.7 eV onto the 
sum of the light neutrino masses |47] . This implies 

m;<0.23eV (33) 

for each mass eigenvalue. Analogous to the limit from Oz//?/? decay in the previous 
section, this bound is modified substantially by the RG evolution. This is shown in 
Fig. El for the eigenvalue m^. As discussed in Sec. I2.3.4| the running of the mass 
eigenvalues is not sensitive to the mixing parameters in the SM, but it depends on the 
Higgs mass. In the MSSM, the variation of the phases causes a slight modification of the 
running, but its order of magnitude is only a few percent even for the large tan j3 used in 
the plot. The influence of #13 is negligible. Interestingly, the evolution of the sum of the 
mass eigenvalues is virtually independent of the mixing parameters for nearly degenerate 
neutrinos both in the SM and in the MSSM. This can be explained by considering the 
sum of the RGEs (JT3J). For mi ~ m 2 ~ m 3 , the terms proportional to add up to 1, 
with small corrections of the order of Am f m and #13. 
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Figure 12: Extrapolation of the upper limit on the neutrino mass from WMAP, rrii < 0.23 eV to higher 
energies, represented by the running of the mass eigenvalue 7713. The SM curves correspond to Higgs 
masses of 114 GeV, 165 GeV and 190 GeV (from bottom to top). In the MSSM, a light Higgs mass of 
120 GeV is used. 

3.3 Constraints on Neutrino Properties from RG Effects 

One may wonder if deviations from 13 = and #23 = vr/4 exist which are the conse- 
quence of radiative corrections. Let us assume therefore that #13 = or #23 = tt/4 are 
given by some high-energy model. Low-energy deviations from the exact values are then 
RG effects, which can be compared to the sensitivities of future experiments. Therefore 
we investigate in a model-independent way the size of RG corrections to 6* 13 and 623 
from the running of the effective neutrino mass operator between the see-saw scale and 
the electroweak scale. 

3.3.1 Corrections to 13 

As pointed out in Sec. I2.3.2[ it is a rather good approximation to assume #13 ~ const, 
in Eq. (jHJ), which leads to an RG evolution with a constant slope depending on the 
Dirac CP phase 5 and the Majorana phases tpx and y?2- Therefore, let us first apply the 
naive estimate (}2*2*|) explicitly to the change of 813 in the MSSM for nearly degenerate 
neutrinos. In this case, the enhancement factor m 2 /Am 2 tm leads to a generic change 
of 813 under the RG that exceeds the detection limit of future experiments even for 
moderate values of tan/3. For example, mi = 0.1 eV and tan/3 = 30 yield a change 
in sin 2 2^x3 of Asin 2 26* 13 ~ 0.5 • 10~ 2 , which is further enhanced by a factor of 4 if the 
Majorana phases are aligned properly. 

In order to obtain a more detailed picture, we now apply Eq. (JHJ) to calculate the 
RG correction to the initial value 813 = between some high energy scale Mi, where 
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neutrino masses are generated, and low energy, i.e. 10 2 GeV. In this case the initial 
value of the Dirac phase 5 is determined by the analytic continuation Eq. ([25)1 . For 
the examples we take M\ = 10 12 GeV. The approximate size of the RG corrections 
to sin 2 2#i3 in the MSSM is shown in Fig. ED In the upper diagram it is plotted as 
a function of tan/3 and the lightest neutrino mass ni\ for constant Majorana phases 
(fx = and (fi>2 — tt- The lower diagram shows the dependence of the corrections on ipi 
and (fi>2 for tan/3 = 50 and mi = 0.08 eV in the case of a normal mass hierarchy. The 
diagrams look rather similar for an inverted hierarchy. Analytically, the pattern of the 
upper plot is easy to understand, and for the lower one there is a simple explanation as 
well. Consider partially or nearly degenerate neutrino masses. Then Eq. (jOJ yields to a 
reasonably good approximation 

6n ~ sin 26>i2 sin 26» 23 -r — — [cos(y?i - 5) - cos(y? 2 - 5)] 

. (fx + - 25 . fix - fi 2 

oc sin sin . d4 

2 2 v ; 

Applying an analogous approximation to Eq. ([25)1 . it can easily be shown that the first 
term in the second line is always ±1, so that the running is completely determined by 
the difference of the Majorana phases. This leads to the diagonal bands in Fig. ED in 
particular the white one corresponding to fii — fi>2 = 0. If one starts with a small but 
non-zero #13, which allows an arbitrary 5, it turns out that the RG evolution quickly 
drives 5 to a value satisfying Eq. (J2SJ), so that the final pattern of Fig. [TH] is unchanged. 

Planned reactor experiments [SU] and next generation superbeam experiments [UT]m^] 
are expected to have an approximate sensitivity on sin 2 2# 13 of 10~ 2 . From Fig.^Jwe find 
that the radiative corrections exceed this value for large regions of the currently allowed 
parameter space, unless there are cancellations due to Majorana phases, i.e. <pi ~ fi>2 
(which might be due to some symmetry). If so, the effects are generically smaller than 
10 -2 as can be seen from the lower diagram. Future upgraded superbeam experiments 
like JHF-HyperKamiokande have the potential to further push the sensitivity to about 
10 -3 and with a neutrino factory even about 10~ 4 might be reached. 

From the theoretical point of view, one would expect that even if some model pre- 
dicted #i3 = at the energy scale of neutrino mass generation, RG effects would at least 
produce a non-zero value of the order shown in Fig. ED Consequently, experiments with 
such a sensitivity have a large discovery potential for 9 13 . We should point out that this 
is a conservative estimate, since if neutrino masses are e.g. determined by GUT scale 
physics, model-dependent radiative corrections in the region between Mi and Mqut con- 
tribute as well [HlinilSnilSlllSSllSSl and there can be additional corrections from physics 
above the GUT scale [HZj- On the other hand, if experiments do not measure #13, this 
will improve the upper bound on #13. Parameter space regions where the corrections are 
larger than this bound will then appear unnatural from the theoretical side. 
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Figure 13: Corrections to #13 from the RG evolution between 10 2 and 10 12 GeV in the MSSM, calculated 
using the analytical approximations with initial conditions #13 = and LMA best-fit values for the 
remaining parameters. The upper diagram shows the dependence on tan (3 and on the mass of the 
lightest neutrino for the case of a normal mass hierarchy and phases ipi = and tf2 = tt. In the lower 
diagram the dependence on the Majorana phases ip\ and ip2 is shown for tan/3 = 50 and mi = 0.08 eV. 
The contour lines are defined as in the upper diagram. In order to apply Eq. JjJJ to the case #13 = 0, 
where 5 is undefined, the analytic continuation of Eq. I|25|) has been used. 
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3.3.2 Corrections to 2 3 

We now consider the RG corrections which induce a deviation of #23 from n/4, even 
if some model predicted this specific value at high energy. We apply the analytical 
formula fllOD with a constant right-hand side in order to calculate the running in the 
MSSM between Mz and the see-saw scale, which we take as M\ = 10 12 GeV for our 
examples. As initial conditions we assume small #13 at Mi and low-energy best-fit 
values for the remaining lepton mixings and the neutrino mass squared differences. In 
leading order in 13 , the evolution is of course independent of the Dirac phase S. 

The size of the RG corrections in the MSSM is shown in Fig. El From the upper 
diagram it can be read off for desired values of tan j3 and the lightest mass eigenvalue mi 
in an example with vanishing Majorana phases. The lower diagram shows its dependence 
on the Majorana phases (pi and (p% for tan/5 = 50, mi = 0.1 eV and a normal mass 
hierarchy. The diagrams look rather similar in the case of an inverted hierarchy. The 
effects of the Majorana phases can easily be understood from Eq. (110(1 . In the region with 
ipi tp 2 ~ vr (again, this might be, e.g., due to some symmetry), both |m 2 e 1I/P2 + m^ 2 
and I mi e l{pl +m 3 | 2 are small for quasi-degenerate neutrinos, which gives the ellipse with 
small radiative corrections in the center of the lower diagram. Such cancellations are 
not possible with hierarchical masses, but the RG effects are generally not very large in 
this CclSC, clS shown by the upper plot. 

Even if a model predicted #23 = 7r/4 at some high energy scale, we would thus expect 
radiative corrections to produce at least a deviation from this value of the size shown in 
Fig. El so that experiments with such a sensitivity are expected to measure a deviation 
of #23 from 7r/4. The sensitivity to sin 2 2^23 of future superbeam experiments like JHF- 
SuperKamiokande is expected to be approximately 1% (see e.g. [HE])- This can now be 
compared with Fig. El We find that the radiative corrections exceed this value for large 
regions of the currently allowed parameter space, where no significant cancellations due 
to Majorana phases occur. This means that tpi and ^2 must not be too close to 71. 
Otherwise, the effects are generically smaller as can be seen from the lower diagram. 
Upgraded superbeam experiments or a neutrino factory might even reach a sensitivity 
of about 0.5%. As argued for the case of 6>i 3 , if experiments measure #23 rather close to 
7r/4, parameter combinations implying larger radiative corrections than the measured 
deviation will appear unnatural from the theoretical point of view. 
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Figure 14: Corrections to 6> 2 3 from the RG evolution between 10 2 GeV and 10 12 GeV in the MSSM, 
calculated from the analytical approximation Eq. (|10f> with initial conditions 623 — tt/4, small #13 = 
and LMA best-fit values for the remaining parameters. The upper diagram shows the dependence on 
tan (J and on the mass mi of the lightest neutrino for the case of a normal mass hierarchy and phases 
<Pi — f2 = 0. In the lower diagram the dependence on the Majorana phases ipi and ipi is shown for 
the example tan/3 = 50 and m\ =0.1. Note that for small $13 the results are independent of the Dirac 
phase to a good approximation. 
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4 Conclusions 



We have derived compact expressions which allow an analytical understanding of the 
running of neutrino masses, leptonic mixing angles and CP phases in the SM and MSSM. 
The results are given directly in terms of these quantities as well as gauge and Yukawa 
couplings, and especially for a small angle 13 the expressions become very simple, 
even when non- vanishing CP phases are present. We have extensively compared those 
formulae to numerical results and we have found that the RG evolution of the physical 
parameters is described qualitatively, and to a reasonable accuracy also quantitatively, 
very well. We have shown that Dirac and Majorana CP phases can have a drastic 
influence on the RG evolution of the mixing parameters. We have reproduced and 
illustrated some effects that were previously described in the literature. As a particularly 
interesting example, we have discussed the radiative generation of the Dirac phase from 
the Majorana phases. Besides, we have derived new results, for example concerning the 
running of the CP phases. Even though the RG effects for the mixing parameters in the 
SM are rather small, the RG effects for the masses are not, and have to be taken into 
account in any careful analysis which relates high and low energy scales. In the MSSM, 
especially for large tan /3, the evolution of the mixings and phases can be large. 

The RG evolution has interesting phenomenological implications. In the case of 
leptogenesis, we have estimated the corrections which arise if the running is appropriately 
taken into account and found that the mass window for baryogenesis is likely to shrink 
when those corrections are considered. In order to simplify the inclusion of RG effects 
in future calculations, we provide the relevant information of the mass parameters at 
the leptogenesis scale. Furthermore, we investigated the extrapolation of the upper 
bounds on the neutrino mass scale from 0uf3f3 decay experiments and WMAP to higher 
energy scales, where they become restrictions for model building. Experimentally one 
finds 9 2 s ~ 7r/4, 13 ~ 0. The deviations from 7r/4 and zero may have a radiative 
origin and we calculated therefore in a model-independent analysis the RG corrections 
to 023 = 7r/4, #i3 = 0. With future precision experiments this may lead to interesting 
insights into model parameters. 

To conclude, we have obtained analytic formulae which are a useful tool to under- 
stand the RG corrections, relevant whenever parameters at two different energy scales 
are compared. This has been demonstrated in the phenomenological applications. 
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Appendix 



A Definition and Extraction of Mixing Parameters 
A.l Standard Parametrization 

In this section we describe our conventions and how mixing angles and phases can be 
extracted from mass matrices. For a general unitary matrix we choose the so-called 
standard-parametrization 

U = diag(e iS %e i5 ^e i5 ^ ■ V • diag(e^ l/2 , e~ iLp2/2 , 1) 



where 



C12C13 



S12C13 



sue 



-id 



V = 



-C23S12 - S23Sl3Cl2e l5 C23C12 - S 2 3Sl 3 Sl2e 1<5 S23C13 



AS 



AS 



(A.l) 



(A.2) 



-S23 s 12 — C23Si3Ci2e — S23 c 12 — C23Si3Si2e C23C13 

with Cij and Sij defined as cos 9ij and sin%, respectively. 
A.2 Extracting Mixing Angles and Phases 

In this standard-parametrization, the mixing angles #13 and #23 can be chosen to lie 
between and |, and by reordering the masses, #12 can be restricted to < 612 < \ . 
For the phases the range between and 2n is required. In order to read off the mixing 
parameters, we use the following procedure: 



1. 6>i3 = arcsin {\U 13 \). 

arctan 



2. 9 12 — 

3- #23 — 

4. 5, 

5. 5 T 



(■ 


u 12 







arctan 





U23 


( 


U33 



if two 

else 

ifC/337^0 
else 



6. S 



arg(C/ 23 ) 
arg(C/ 33 ) 

/ u*u u u 

— arg 
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Cl2 C 13 C23 S13 



+ c 



12 c 23 s 13 



Sl2 S23 

where i, j G {1, 2, 3} and i 7^ j. 



V 



30 



7. 5 e = arg(e i5 f/ 13 ) 

8. ^ 1 = 2arg(e i5e f/* 1 ) 

9. ip 2 = 2^g(e^U* 12 ) 
Here we used the relation 

UijUijUjiUjj = C12 Ci 3 C 23 Si3 (e _1<5 Si 2 S 23 - C12 C 23 S13) , 

which holds for i,j e {1, 2, 3} and z ^ j. Note that this relation is often used in order 
to introduce the Jarlskog invariants [HI] 

J CP = i|M^W2iC^)| = i|Im(^i3^3i£4)l 

= j I Im ( ^22 ^2 3 C/32 C^33 ) I = g l C 12 C 13 C 23 SU1 5 S 12 S 13 S 23 | • (A.3) 

For the sake of a better numerical stability, one can choose any of the three combinations. 
In particular, if the modulus of one of the is very small, it turns out to be more 
accurate to choose a combination in which this specific Uij does not appear. 

A.3 Leptonic Mixing Matrix 

Since the effective neutrino mass matrix is symmetric, it can be diagonalized by a unitary 
matrix U v , 

Ul m u U u = diag(mi, m 2 ,m 3 ) . (A.4) 

The form of U depends on a prescription how to order the mass eigenvalues. In order 
to obtain a mixing matrix which can be compared with the experimental data, the 
choice of the prescription is somewhat subtle. From experiment we know that there 
is a small mass difference, called A?7ig 0l = raf — m|, and a larger one, referred to as 
Am^ tm = m\ — mj. By convention, the masses are labeled such that i,j ^ 3 while 
either k or £ equals 3. The different schemes are depicted in Fig. ED The mass label 2 
is attached to the eigenvector with the lower modulus of the first component. We are 
doing this since we want to read off a mixing angle #i 2 less then 45°. 

The neutrino mixing matrix Umns can then be read off in the following way: 

1. Diagonalize Y}Y e by U e , i.e. Y e -> U\ ■ Y% -Y e -U e = diag (y%, y*, yfj where yj are 
positive for / e {e, (/,, r}. 

2. Change the basis according to m u — > m' u = f/J • m v ■ U e . 

3. Diagonalize m' u : m' u — > ?7^ NS • m' v ■ ?7mns — diag(mi, m 2 , m 3 ) where m, > 0. 

Then ?7mns contains the leptonic mixing angles which can be read off as described in 
Sec. IA.2I Note that mi < m 2 < m 3 is not necessarily fulfilled, as we already mentioned 
before (cf. Fig. fTHjl . 
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■m 2 /i 

■ m 1 12 m 3 



(a) Normal mass hierarchy. (b) Inverted mass hierarchy. 

Figure 15: The normal and inverted mass hierarchy. 

B Derivation of the Analytical Formulae 

To derive the RGEs for the mixing parameters, we follow in general the methods of [70J. 
The RGE for k reads 

16vr 2 ^ = aK + P T K + KP, (B.5) 

where all terms with trivial flavour structure are absorbed in a. k can be diagonalized 
(in the basis where Y e is diagonal) by a unitary transformation, 

U{t) T K (t)U(t) = D(t) = — di&g (m 1 (t),m 2 (t),m 3 (t)) . (B.6) 



We hence obtain 



4- (U* D U ] ) = U* DU ] + U* DU ] + U* DU ] 
at 

™ (aU* DU ] + P T U* DU ] + U* DU ] P) . (B.7) 
IQtt 2 x ' 

Multiplying with U T from the left and with U from the right yields 

U T U* D + DU ] U + D = \aD + P' T D + D P'] , (B.8) 

where we have introduced P' = U< P U. The next step is defining an anti-Hermitian 
matrix T by 

^-U = UT . (B.9) 

dt v ; 
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With this definition, we find 

t> = — — (aD + P' T D + D P') — T* D + DT , (B.10) 

where the anti-hermiticity of T was used. Since the left-hand side of this equation is 
diagonal and real per definition, the right-hand side has to possess these properties as 
well, 

™* = 77T~2 ( am i + 2 P'u m i) + (Tu — T-*) rrii . (B.ll) 

l07T^ 

Note that here and in the following equations, no sum over repeated indices is implied. 
The second bracket is purely imaginary, hence it has to cancel with the imaginary part 
of the first one, 

2ImT« = -^(Ima + 2Im^) , (B.12) 

and we further confirm eq. (15) of [11J, which translates with our conventions to 

16vr 2 rhi = (Re a + 2 Re P(J m, . (B. 13) 

Eq. dBHH) differs from Eq. (19) of [TT], where the imaginary part of a is not present; 
however, this difference is irrelevant in the SM and the MSSM, where a is real. By 
comparing the off-diagonal parts of (jB.lOjl we find 

m, Tij - 7* m, = -^ [P[J m, + m % PQ . (B. 14) 

Adding and subtracting this equation and its complex conjugate, we obtain for % 7^ j 

„ m j ReP' { + rrii ReP/, 

16tt 2 ReT i7 - = J - ^ , B.15a 

rrii — trij 

„ rrij ImP', + rrii ImP-, 
16tt 2 ImT- = 3 - ^ . B.15b 

rrii + mj 

Let us now focus on Hermitian P, which implies Hermitian P', for a moment. Using 
RePjj = Re Pi* = ReP/^ and an analogous relation for ImP/-, we obtain in this case 

16vr 2 ImT i7 - = _ mi ~ m i Im P' (B.16a) 

3 rrii + mj 13 v ' 

lG^ReT,, = - Wt +mi ReP' (B.16b) 
rrii — m j 

In order to obtain the renormalization group equations for the mixing angles, we use 

una, 

U ] U = T. (B.17) 
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Inserting the standard parametrization (jA.ljl . we can express the left-hand side of (jB.17|) 
in terms of the mixing parameters and their derivatives. Now we can solve for the 
derivatives of the mixing parameters. Note that due to the separation of the evolution 
of the mass eigenvalues in equation (|B.13J) . we have reduced the number of parameters 
from 12 to 9. The discussion so far has been very similar to the one of [H]. There, 
the RG evolution of the mixing parameters is expressed in terms of the mixing matrix 
elements and P' . 

In order to obtain rather short and more explicit formulae, which are e.g. useful for 
deriving the approximations of Sec. 12.11 we now consider (|B.17|) and label the mixing 

parameters as 

{&} = {012, 013, 023, 5, S e , 8ft, 5 T , <fi, ip 2 } . (B.18) 

We observe that the left-hand side of ()B.17|) is linear in Therefore, by solving the 
corresponding system of linear equations, we can express the derivatives of the mixing 
parameters by the mixing parameters, the mass eigenvalues and the Yukawa couplings. 
The resulting formulae are still too long to be presented here but can be obtained from 
the web page http://www.ph.tum.de/~mratz/AnalyticFormulae/. 

Finally, let us record that only the moduli of Uij enter into the diagonal elements of 
P', if P is diagonal, P = diag(Pi, P2, P3) (which is the case in the SM and MSSM in the 
basis we have used in the main part), since 

p li = ^(j^hPjkUki = Yl u *ji P 3 5 ohUki = \ U 3i\ 2p i ■ (B.19) 

jk jk j 

Consequently, the evolution of the mass eigenvalues does not directly depend on the 
Majorana phases, as claimed in Sec. 12.3.41 
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